Temperature dependent shifts in collective excitations of a gaseous Bose-Einstein 

condensate 
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The moment method is applied to the quantum theory for a trapped dilute gas, obtaining equa- 
tions for the evolution of the cloud. These equations proof the existence of undamped oscillations in 
a two-dimensional harmonic trap with radial symmetry. In all other cases we can, with physically 
reasonable approximations, solve the equations and find the frequencies of the collective modes for 
all temperatures, from the pure Bose-Einstein condensate up to the normal gas. 

PACS numbers: 



Ever since the first achievement of a weakly interacting 
Bose-Einstein condensate in laboratory [El 0, 0], one of 
the basic tools for its study were the oscillation modes of 
the condensed cloud [H |j. The measure of the frequen- 
cies of these oscillations at very low temperatures, and 
the comparison with several theoretical predictions |], |7j 
were the first important backing to the Gross-Pitaevskii 
equation, the mean-field model which describes the con- 
densed cloud in the zero temperature limit. 

Since those first works, there has been a growing inter- 
est on the study of the intermediate regimes, which are 
far from the zero temperature limit and closer to the crit- 
ical temperature T c at which condensation begins. For 
T > 0.5T C the collective modes of the three-dimensional 
condensate suffer both a frequency shift and a strong 
damping Q|. These deviations from the mean field the- 
ory are due to the uncondensed or normal component 
of the gas. The role of the the thermal cloud has been 
partially explained using self-consistent solutions of the 
many-body theory for the condensate ||, [|, theory of 
relaxation phenomena, linear response theory, etc [See 
Ref. jlO) for a list of references] . Nevertheless the theo- 
retical and experimental studies of these excitations and 
of their damping are far from complete, and although 
the sum of all works gives some insight into the behavior 
of the dilute gas, there is no macroscopic theory which 
systematically predicts both the finite temperature shifts 
and the damping of the excitations. 

In this work we attempt a different approach to the 
study of collective modes. Our goal is not to consider 
the condensed and normal components separately, but 
to treat the system as a whole. We will start from first 
principles, studying the quantum Hamiltonian of the di- 
lute gas and deriving equations for a few relevant ob- 
servables: the widths of the cloud, the kinetic energy on 
each direction, etc. We will prove that these equations 
close exactly in the two-dimensional case with radially 
symmetric trap. This exact closing reveals the existence 
of undamped radial excitations of the cloud, with a uni- 
versal frequency u>2D = 2w which is twice the frequency 
of the trap. In the last part of this work we will return 



to the moment equations for arbitrary geometry and we 
will close them using a reasonable approximation for the 
interaction energy as a function of the widths. With this 
hypothesis the dynamical equations can be solved, giving 
us the temperature dependence of the excitation frequen- 
cies in the cloud as a whole. Finally this dependence is 
compared to other works || and to experiments p|. 

The model.- The theory for a dilute gas of bosons in 
a harmonic confinement is developed up from a simple 
Hamiltonian 



H = / *(x) 



--A + v(x,;) 



J ^¥(x) + ¥(x)t*(x)¥(x)d B a;. 



(1) 



The trapping potential V(x, t) — \ J2i &i(t) 2x i > ma Y 
have any symmetry, and it may be even subject to 
time dependent perturbations. The coupling constant 
U = 4irh 2 a/m measures the interaction among bosons in 
terms of the s-wave scattering length a and the contact 
interaction Vbosons ( x , y) = a5(x— y). Finally, the Hamil- 
tonian has an important hidden parameter which is the 
dimensionality of the space, n € {1, 2, 3}, and that has a 
crucial role on the behavior of the cloud. 

The Hamiltonian H must be understood as an opera- 
tor on a Fock space, where the creation and destruction 
operators obey the commutation rules 

[*(x),*t(y)] = <5(x-y), [*(x),*(y)]=0. (2) 

Below a critical temperature, T c , the dilute gas experi- 
ences a phase transition which leads to a macroscopic 
population of the ground state of Ti.. By neglecting the 
number of atoms in the thermal cloud one may approx- 
imate *(x) ~ a//V0(x) + 5\&(x), where N is the total 
number of particles of the cloud. The c-number 0(x) 
constitutes the order parameter of a mean-field theory, 
and it obeys the so-called Gross-Pitaevskii equation 
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V(x,t) + UN\(f>\ 2 



(x,t). (3) 
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FIG. 1: Frequencies for monopolar (i>_ ,+), and quadrupolar 
(1/2) excitations of a dilute gas, versus normalized tempera- 
ture t = T/T c . We plot the theoretical predictions for (a) 
spherical trap Ui = 1 and (b) pancake trap 7 = 8 interpo- 
lating the interaction energy between the Thomas-Fermi and 
ideal limits. All figures have been adimensionalized. 



In the Thomas-Fermi limit for a stationary trap, the ki- 
netic energy of the wavefunction is neglected with respect 
to the interaction energy. Then the system admits sta- 
tionary solutions |<?!>t.f| 2 oc (fi — V(x))/UN which allows 
for a self-consistent determination of the chemical poten- 
tial in terms of the trapping strength and of the number 
of particles, N. Using (\>tf one may study the monopo- 
lar and quadrupolar excitations of the condensate fl and 
compare them to the experimental results 0). The con- 
clusion is that Eq. (|J) provides a good description of the 
condensate when the thermal cloud is undetectable, while 
for higher temperatures the collective modes experience 
a strong, temperature dependent shift B. This shift has 
been sometimes explained by changing the effective num- 
ber of particles in the condensate, N, and recomputing 
the normal modes, or by performing a numerical study 
of the Hamiltonian in various self-consistent frameworks 

The moment equations.- The word "moment" in 
this context refers to the expected value of an observable. 
If we work with the mean-field theory these moments are 
calculated by averaging over the order parameter. For 
instance the center of mass of the condensed cloud is 



X = 



(x)x(f>(x)d n x, 



(4) 



and it follows an equation which is exact, Xi = 
—uji(t) 2 Xi. These equations and similar ones for other 
expected values are powerful tools both in the study of 
Bose-Einstein condensates and in the study of nonlinear 
Schrodinger equations in general jll], [l2| . 

We will apply an equivalent technique to the study of 
the quantum Hamiltonian (^). First of all we adimen- 
sionalize 7i, using as fundamental units the transverse 



frequency of the trap, uj± — ^u x ujy, and the radial size 
of this trap, do = H/ ' ■ s fnuJT\_. This way, the Hamiltonian 
becomes, up to a global factor, 



H = J *(x)t [-±A + Er hui{t?x1] *(x)d" 



I *(x) t * (x) f * (x) 4? (x) d n x, 



(5) 



where u>i = aji/uj± and g = Ana/aQ. Next we define a set 
of one-particle operators 



X 
P 

Wi 
Bi 
Ki 
J 



W(x)x9(x)d n x, 
¥ t (x)(-tV)*(x)d n ar, 
¥(x)x 2 9{x)d n x, 
¥(x){~ixA)9(x)d n x, 
-\ J &(x)df*(x)d n x, 

■* t (x)* t (x)*(x)*(x)d"x, 



(6a) 
(6b) 
(6c) 
(6d) 
(6e) 
(6f) 



I J [{d^)H 2 - {^f(d^f]d n x. (6g) 



In these equations di represents a derivative with respect 
to Xi. From top to bottom, we have defined the center 
of mass, its moment, the width of the cloud, the rate of 
change of the width, the kinetic energy along the i-th 
direction, the interaction energy and the rate of change 
of J along the i-th axis. 

We use the Heisenberg picture to study the evolution 
of these operators. In this image the density matrix of 
the system dictates some initial conditions, p, and the 
time evolution is carried by the observables. Thus, any 
observable A without explicit time dependence follows 
the simple equation, i-^A = [A, H\, and its expected 
value is simply A = Ti{A(t)p}. 

The first and most important equation is that of the 
destruction or creation operators, 



i&tf (x) = [-|A + !wi(t) 2 a£ + U^y] 9. (7) 

Combining Eqs. (||)-(0), and using the fact that our oper- 
ators lack an explicit time dependence, we find the equa- 
tions for the center of mass 



rf 2 X 
dt 2 



plus a new set of equations for the widths 
D t Wi = B h 

D t Bi = 4K l -2ij l (t) 2 W l + 2J, 
D t Ki = -^{tfBi-Fi, 

DtJ = Ei-P*. 



(8) 



(9a) 
(9b) 
(9c) 

(9d) 
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where D t denotes derivative with respect to time. 

The first important remark is that Eqs. are all 

exact and describe the evolution of one-particle operators. 
Nevertheless, one may now take expected values around 
all terms in Eqs. (||)-(^|) and the same equations will hold 
for the expected values 



X = (X) = Tr{Xp}, 



(10) 



which is on what we focus from now on. The second 
remark is that there are less equations than moments, 
which seems to prevent us from completely determining 
the behavior of the cloud for arbitrary conditions. 

Undamped oscillations.- There is an important sit- 
uation in which the moment equations may be solved, 
and that is a two-dimensional condensate in a radially 
symmetric trap, {n = 2,Wi = oj 2 = oj}- In this geometry 
we define a total width, r = \/Wi + W2, and obtain the 
following equation 



d 2 r , N9 2M 

This equation includes a conserved quantity 



(11) 



M = {K 1 +K 2 + J)(Wi + W 2 ) - \{B X + B 2 ) 2 , (12) 

that must be determined up from the initial conditions. 

In the case of an static trap, the equation for the width 
has an equilibrium point given by R = 2M/u> 2 . The 
oscillations around this point have a proper frequency, 
ll>2d = 2w, which is independent of the shape of the 
cloud. Therefore, the cloud may be at any temperature 
and bear excitations of any multipolarity, but the total 
width of the cloud will always oscillate with the universal 
frequency L1J2D of the monopole mode. 

Furthermore, since Eq. |ll]) is exact and it is con- 
servative, any perturbation of the width of the two- 
dimensional condensate persist eternally, and the damp- 
ing of these oscillations can only be explained in terms of 
higher order contributions to H , such as three- and four- 
body collisions, interaction with the environment and an- 
harmonicities of the trap. 

One may compare Eq. (|ll]) with a similar equation 
arising from scaling arguments pl|. Those arguments 
are based on a symmetry of Eq. (0) such that it admits 
infinite many rescaled solutions 



*(x,t) = i* (x/A,t)e- ,;a;2A /( 2A \ 
A 



A = -cu(t) 2 X + X 



-3 



= A 



(13a) 
(13b) 



for any given solution x I / o( x j t). The parameter A can thus 
be related to the width of the cloud, but the problem with 
scalings is that they implicitly assume a certain type of 
evolution for the operator, while Eq. ( |TT| ) makes no such 
ansatz and it includes further information about final 
shape of the collective mode, thanks to M. 



Asymmetric traps.- We now want to study other 
dimensionalities and other geometries of the trap. An 
inherent limitation in the study of the moments is that 
these hierarchies rarely give us closed equations. One 
must, at some point make a reasonable approximation 
to estimate some of the unknowns and get a tractable 
model. Roughly, we will assume that the self-interaction 
energy is inversely proportional to the effective volume of 
our n-dimensional cloud, V = (HiWi) 1 ' 2 , with a constant 
Jo that depends on the initial data: 



J ~ J (n,^)" 1/2 , F t ~ -\Bt J/Wi. 



(14) 



With Eq. (|l4| ) we can close the moments hierarchy, 
obtaining 2n differential equations for 2n moments 



d Xi 

dMi 
dt 



J- 



-LUi(t) 2 Xi 

M 
dt ' 



2M, 



J_ 



(15a) 
(15b) 



with three new variables Mi — KiWi — \B 2 . 

The motivation for the ansatz (|14j) is multiple. First 
it arises from scaling considerations. If we assume that 
{Wi . . . W n } are the only independent variables of our 
model, and we rescale our model according to Xi — > AjXj, 
the destruction operator, the widths and the interaction 
change as * — > ^ /liiXi, Wi -> X 2 W l , and J -> J/IFA;, 
which immediately leads to Eq. (14). 

The other motivation for Eq. ( |l4|) is that it works in 
the mean field theory too. Making simulations of Eq. (j^) 
with wavefunctions of almost any reasonable shape, it is 
readily seen that Jo is conserved up to a 10%, and even 
in the worst cases the model equations (|l5|a)-(|l5|b) are 
extremely accurate [Q. 

Collective modes.- Eqs. (|l5|a)-(|l5|b) model the evo- 
lution of the widths of a dilute gas at any temperature. 
They do not make any distinction between condensed 
and normal components and the role of temperature is 
to change the equilibrium values of M, and J. In order 
to study the collective modes of the condensate it is eas- 
ier to go back to Eqs. (^a)-(gd) and linearize around the 
equilibrium points {Wi, Bi,Ki,J}. To first order in the 
small variables {wi,b i} fcj}, with the ansatz (ph, we get 



Wi 



Ah - 2u 2 Wl -Y.W- 



j 3 



1 2, J . 

2 1 1 2Wi 



(16a) 
(16b) 

(16c) 



There are two important limits in these equations. In 
the ideal gas limit the interaction energy can be ne- 
glected, J = 0. This means that each width decouples 
from the rest, and that there are three normal modes 
which oscillate with frequencies Vi = 2u>i. In the the 
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Thomas-Fermi limit we rather assume K j <C J and just 
keep the interaction energy and the trap. Then the 
equilibrium point is given by J = ujfWi. For a three- 
dimensional trap with axial symmetry, uj\ = ll>2 = F 
^3 = 7; the excitation frequencies become 



^ 7 + 2±iV9r 



I67+I6, (17) 



in full agreement with Q. The first frequency corre- 
sponds to a m = 2 mode in which the w x and w y widths 
oscillate with opposite phases. The second two frequen- 
cies correspond to two m — modes such that the trans- 
verse shape of the condensate is preserved and w x = w y . 
There are also three Goldstonc modes, v — 0, which arise 
from the scaling symmetry of Eq. (|^). 

Temperature dependent frequencies.- If we now 
focus on the temperature, the ideal gas limit and the 
Thomas-Fermi limit qualitatively describe the T > T c 
and T = limits, respectively. In between these ex- 
tremes the amount of condensed and normal clouds vary, 
inducing changes on the values of {Wi, J} which cause a 
continuous evolution of the oscillation frequencies, from 



{2cji} to values which are close to Eq. ( 

In the radially symmetric case, Eq 
different oscillation frequencies 



17). 

(16) gives us two 



z/_, 2 = w\/4 - 2P, y + = wV4+2P, 



(18) 



where P = J/(ujR 2 ) depends on the actual shape of the 
cloud. A crude estimate for the interaction energy |l^] 
assigns P = 1 and P = for the T = and T = T c limits. 
In between these limits the interaction energy is supposed 
to scale as P = {N /N) 2 / 5 = (1 - t 3 ) 2 / 5 , where N is the 
number of particles in the condensate, and t = T/T c is 
the temperature relative to the phase transition due to 
the Bose- Einstein condensation [|l0| . This oversimplified 
model provides the same frequencies as a condensed cloud 
with No particles and no thermal cloud, and it agrees well 
with self consistent calculations in the Popov theory || . 

As an example, in Figure |l|(a) we plot the oscillation 
frequencies as a function of the relative temperature, t, 
for a spherically symmetric trap. Better estimates of 
the excitation frequencies require better estimates of the 
functions {J(T),Wi(T)}. That task is significantly sim- 
pler than full calculations of the excitation spectrum for 
N particles, but it is yet an open problem. 

Finally we have compared our predictions with the ex- 
perimental results of D. S. Jin et al. The only dif- 
ference is on the quadrupolar excitations, whose 
experimental frequency seems to decrease when the tem- 
perature increases. This point was already notice in ||] 
but now it becomes evident that the deviation could be 
due to studying the width of the condensate and of the 
thermal cloud separately. We believe that this appar- 
ent contradiction deserves further experimental investi- 
gation, preferably by studying the oscillations in clouds 
with more atoms and various geometries. 



Summary and discussions.- We have studied the 
collective modes of a trapped dilute gas. By computing 
the evolution of several operators, we have obtained dif- 
ferential equations for the center of mass and the widths 
of the cloud. In the two-dimensional case with radially 
symmetric trap we have predicted the existence of un- 
damped oscillations of the radial width with frequency 
u>2D = 2w twice the frequency of the trap. This predic- 
tion can be confirmed in current experiments |l4| . 

For other symmetries and dimensionalities we recov- 
ered both the Thomas- Fermi limit and the ideal gas limit. 
We also combined a thermodynamical estimate of the in- 
teraction energy with our dynamical equations to obtain 
the temperature dependence of the collective modes in 
various traps, with similar results as previous numerical 
works Q and experiments Q . 

This work offers new techniques that can be applied 
to the study of trapped dilute gases in general and to 
Bose-Einstein condensation in particular. We suggested 
a new way to perform experiments, by studying the gas 
cloud as a whole, and we also offered a simple way to 
relate the normal modes of the cloud to the interaction 
energy of the gas, without the need to perform difficult 
self-consistent calculations. Finally, with our analysis it 
becomes clear why simple models such as scalings |l^] 
have been so useful in analyzing current experiments. 



[7; 
[§ 
[9 
[10 

[11 

[12. 

[13 
[14 



M. H. Anderson et al, Science 269, 198 (1995). 

K. B. Davis et al, Phys. Rev. Lett. 75, 3969 (1995). 

C. C. Bradley et al, Phys. Rev. Lett. 75, 1687 (1995). 

D. S. Jin, M. R. Matthews, J. R. Ensher, C. E. Wieman, 
and E. A. Cornell, Phys. Rev. Lett. 78, 764 (1997). 
M.-O. Mewes et al, Phys. Rev. Lett. 77, 988 (1996). 

S. Stringari, Phys. Rev. Let. 77, 2360 (1996); F. Dalfovo, 
L. Pitaevskii, and S. Stringari, Phys. Rev. A 54, 4213 
(1996). 

V. M. Perez-Garcfa et al, Phys. Rev. Lett. 77, 5320 
(1996). 

D. A. W. Hutchinson, E. Zaremba, and A. Griffin, Phys. 
Rev. Lett. 78, 1842 (1997). 

R. J. Dodd, M. Edwards, C. W. Clark, and K. Burnett, 
Phys. Rev. A 57, R32 (1998). 

F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999). 

J. J. Garcfa-Ripoll, V. M. Perez-Garcfa, and P. Torres, 
Phys. Rev. Lett. 83, 1715 (1999). 

J. J. Garcfa-Ripoll, Ph.D. Thesis, Universidad Com- 
plutense de Madrid (2001). 

Y. Kagan, E. L. Surkov, and G. V. Shlyapnikov, Phys. 
Rev. A 55, R18 (1997). 

A. Gorlitz, J.M. Vogels, A.E. Leanhardt, C. Raman, T.L. 
Gustavson, J.R. Abo-Shaeer, A. P. Chikkatur, S. Gupta, 
S. Inouye, T.P. Rosenband , D.E. Pritchard, W. Ketterle, 



arXiv: cond-mat/0104549 



